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The proper quantum plasma treatment ol the electron gas in degenerate stars such as white dwarfs 
provides an additional quantum contribution to the electron pressure. The additional pressure term 
modifies the equation for hydrostatic equilibrium, resulting in the quantum modified Lane-Emden 
equation for polytropic equation of states. The additional pressure term also modifies the expression 
for the limiting Chandrasekhar mass of white dwarfs. An approximate solution is derived of the 
quantum modified Lane-Emden equation for general polytropic indices, and it is demonstrated that 
the quantum corrections reduce the standard Chandrasekhar mass and enhance the white dwarf 
radius by negligibly small values only. 



I. INTRODUCTION 

From a plasma point of view the electron gas of den- 
sity n in degenerate stars such as white dwarfs is a nearly 
collisionless quantum plasma 9 because its central temper- 
ature T is lower than the Fermi temperature Tp, and the 
quantum coupling parameter qq = E c /Ep of the ratio of 
the Coulomb interaction energy E c ~ e 2 ™ 1 / 3 to the Fermi 
energy Ep is much smaller than unity. The white dwarf 
Sirius B has an average mass density of 3 • 10 6 g cm -3 
and a central temperature of T c — 7.6 ■ 10 7 K yielding for 
the temperature ratio 



X^ = (^f 3 f^ =3.65 (nX%r = 55.7 

(1) 

The ratio ((lj indicates that the thermal de Broglie wave- 
length of individual plasma particles 



Here we demonstrate that the proper quantum treat- 
ment of the electron gas based on the Wigner— dis- 
tribution function changes the hydrostatic equilibrium 
equation in white dwarfs leading, however, to a negligi- 
ble modification of the maximum Chandrasekhar mass 
of such systems. In the hydrodynamical equations the 
quantum effect 8 yields an additional pressure term Pq 
to the classical pressure P c , 



P = P c + Pq=P c 



2m 



(Vn^f-n^vV/ 2 ! (3) 



in terms of the electron density n. The additional pres- 
sure term Pq is not listed in Salpeter'&ii account of the 
pressure contributions to a zero-temperature degenerate 
Fermi gas of non-interacting electrons. Here the modi- 
fication to the standard Lane-Emden equation for poly- 
tropic gases and to the limiting Chandrasekhar mass are 
calculated. 
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3.9n 



-1/3 
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is of the same order of magnitude as the mean distance 
between electrons nT 1 ^ . The thermal de Broglie wave- 
length roughly represents the spatial extension of a parti- 
cle's wave function due to the quantum uncertainty. For 
Xb comparable to the interparticle distance, because of 
overlapping electron wave function extensions, individual 
electrons cannot be treated as pointlike particles as in the 
classical plasma description, and quantum effects become 
important. The smallness of the quantum coupling pa- 
rameter <7q assures that collective mean field-effects dom- 
inate over binary collisions because the typical electron 
Fermi energy is much larger than the Coulomb interac- 
tion energy with neighbouring electrons. 
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II. QUANTUM MODIFIED EQUATIONS FOR 
HYDROSTATIC STELLAR EQUILIBRIA 

For a plasma at rest (u = 0) the hydrodynamical equa- 
tion (|A10[) becomes 



= g(x) + -VP C 
P 



n 2 z / vy/ 2 

2mm,i \ p 1 / 2 



(4) 



with the mass density p = nrrii/ Z where mi denotes the 
ion mass. Choosing spherical coordinates we obtain the 
quantum modified equation for hydrostatic equilibria in 
stars as 



= g(r) + 



1 dP r 



p(r) dr 
with the operator 

V 2 = r 



h 2 Z d f Vlp^jr) 
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(5) 
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and the usual gravitational acceleration 



where the dimensionless parameter 



g(r) = Gp(r)/r 2 
due to the enclosed mass 



dfx 
dr 



Anr 2 p{r) 



The mass inside the radius r is given by 



M(r) = / dr ^fj — Att I dr (r') 2 p(r 
Jo dr J 

so that the total mass is 



(7) 
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(9) 



h 2 z 
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characterizes the strength of the quantum plasma modi- 
fications. The term on the right hand side of the QMLE 
equation (|14[) reflects the new quantum contribution from 
the Bohm potential. For i] a = the QMLE-equation re- 
duces to the standard Lane-Emden equation of stellar 
structure. 

The QMLE equation (TT4"]) has been derived for flat 
Newtonian space-time ignoring general relativistic ef- 
fects. For completeness in Appendix B we derive the 
relativistic generalization of the QMLE-equation. 



M = 47r / dr r l p(r) 
'o 



(10) 



Inserting Eq. (J7J into Eq. ((SJ), and differentiating with 
respect to r immediately yields the quantum modified 
stellar structure equation 



d_ 

dr 



1 dP c h 2 Z d /Vjy/3( r ) 



p(r) dr 2mrrii dr \ p 1 / 2 (r) 
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A. Quantum modified Lane-Emden equation 

Following standard procedures^ we adopt the poly- 
tropic equation of state 



P c =Kp 



(l+a)/a 



(12) 



and substitute p = ptfJa (with central density p c ) and 
r = Ax with the constant 



B. Quantum modified Chandrasekhar mass 

In terms of the substituted variables x and y a (x) the 
mass © inside the normalised radius x reads 



Ax 



M(x) = 4:71 I dr (r') 2 p(r ) 
'o 



Jo 

Inserting the QMLE equation (|14l) yields 



M{x) = -M 



> dy a 
dx 



' dx 
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dx 
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with 



(16) 



(17) 
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(13) 



Eq. (jlip then becomes the quantum modified Lane- 
Emden (QMLE) equation for a polytrope of index a 



M = 4tt Pc A 3 



(18) 



The first zero x\ of the solution y a {xi) = of the QMLE 
equation defines the size of the star. In terms of Mo, 
with Z — 2 and mj = Z ■ 1836m, the parameter r\ a can 
be expressed in cgs-units as 



}_d_ ( 2 dy a 
x 2 dx V dx ] + Va 
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Scaling M = 3xM sun = 6x • 10 33 g in solar masses and 
A = ag ■ 10 9 cm with the typical white dwarf mass and 
radius, we obtain very small values of the parameter 
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r, a = 9.2 • 10- J >9X) 



(20) 



The small value of the parameter r\ a reflects the the small 
ratio of the Fermi pressure and Bohm pressure terms (sec- 
ond and third term in Eq. ([5])). By replacing (d/dr) by 
1/L, where L denotes a typical length such as the radius 
of the star, and using the relativistic expression for the 
Fermi pressure P c = (37r 2 ) 1 / 3 ?icn 4 / 3 /4, we find for the 
ratio of the two terms in terms of the Compton electron 
wavelength A c = H/(mc) — 3.86 • 10~ n cm 



T-2 



Xr 



2mP c L 2 (3TT 2 ) 1 / 3 L 2 nV 3 



(21) 



For a white dwarf, typically we have n ~ 10 30 cm 3 , so 
that with L ~ 10 9 cm we estimate T 3 /T 2 ~ 2.5 • 1CT 39 , 
in agreement with the estimate ([211)) . 

While these rough order of magnitude estimates in- 
dicate that the Bohm pressure term is likely a small 
contributor to modifications of the Chandrasekhar white 
dwarf mass limit, nevertheless it is appropriate to figure 
out the effect quantitatively for two main reasons. First, 
it can happen that the rough estimates given above do 
not uncover the total effect when it is worked out quan- 
titatively and so it could happen that a small correction 
has a profound effect on the solution to the non-linear 
equation. Such effects are not unusual in physics and 
one truly needs to work out the details in order to be 
sure that one has not overlooked some subtle component. 
Second, even when detailed calculations show that the 
Bohm pressure effect is small, in accord with the rough 
estimates made, it is then satisfying to note the confirma- 
tion of the rough estimates by detailed evaluation. For 
these two reasons alone it is more than appropriate to 
evaluate the Bohm pressure effect although, as we will 
show, the contribution is indeed small. The procedure 
for evaluating the contribution to the white dwarf mass 
may also be of relevance to other astrophysical problems 
and so, as a basic technique, it is also appropriate to 
spell out the details of how one addresses such non-linear 
corrections. 

Since the pioneering work of Chandrasekhar 1 - 3 sum- 
marized in^ it is well known that the interior of white 
dwarfs is supported by the pressure of completely de- 
generate electrons. With increasing density the pressure 
becomes so high that the degeneracy becomes relativistic 
so that P c = 1.244 • 10 15 (/o//i e ) 4/3 which corresponds to 
a polytrope index of a = 3. In this case Eq. (TT71) yields 
the quantum-modified Chandrasekhar mass limit 
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with 



% = 4.28 • 10 



K 2 



(23) 



The first term in Eq. ((22]) is the standard expression of 
the Chandrasekhar mass limit whereas the second term 
represents the modification from the quantum plasma 
Bohm potential. In order to investigate its effect in en- 
hancing or reducing the standard mass limit one has to 
compute the first zero x\ and the shape of the solution 
2/3 of the QMLE equation. 



III. APPROXIMATE SOLUTION OF THE QMLE 
EQUATION FOR GENERAL POLYTROPIC 
INDICES a FOR SMALL VALUES OF r/ a 

Substituting t = 1/x the QMLE equation (fT4|) reads 



dt 2 



flat 4 d 2 yg /2 
J*l 2 dt 2 

yet 



(24) 



Dropping the index a, i.e. y(t) = y a {t) and rj a = rj, and 
using u(t) = [y{t)] a/2 the QMLE equation reads 



dt 2 



[y- 



rjt 4 d 2 u-, 
~~di 2 ^ 



t 4 



= 



(25) 



A. First-order expansion in r\ 



For very small values of the quantum parameter r\ <C 1 
set 

u(t) = u c (t) + r)Su (26) 
where u c fulfills the standard Lane-Emden equation 



d 2 ul /a -' 2 



^ = 



dt 2 t 4 
The ansatz (|2"6"|) implies 

y{t) = Vc(t) + vhj = (u c + -qSuf^ 



(27) 



M = M(xi,a = 3) = -- 



.1/2 



3/2 



+ ?H u ( 2 -<*y*su 



(28) 
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so that to first order in rj 



Su 



,(a-2)/a 



Sy 



(29) 



Using Eqs. (gU) and (J35J) we find for the QMLE equation 
(1231) to first order in 77 

dfi ^ + V5V -~ c ^ + ? = ° (30) 

Using Eqs. (HU and (22) we derive 



CP r 1 <!■(/,■■- _ 



dt 2 



dt 2 J t 4 



dt 2 J t 4 



+ S«c (Q " 1)/a ^ = (31) 



which, for a known standard solution u c , is a linear 
second-order differential equation for the deviation 6y. 

Both u(t) and the solution u c (f) of the standard Lane- 
Emden equation fulfil the boundary conditions 



3 1 dy c d 



5y 



u c dt 2 

u c (q)Su(q) 
9 4 



(36) 



where we inserted equation (|34p . Eq. (|36l) represents the 
first-order correction in 77 of the Chandrasekhar mass. 

If t c denotes the largest zero of u c (t c ) =0we set t\ = 
t c + St. The boundary condition 

= u(ti) = u c {ti) + rjSu(ti) 
together with the Taylor expansion 



u c (ti) ~ u c (t c ) + (^)tjt = (^)tjt 
at at 



then yields 



so that 



ou(fi) <Sit(t c ) 

or = — 77 — 3 ~ 

' / du c \ ' / cm c \ 

I dt /*c I dt 7*c 



(37) 



implying 



u c (f = 00) = 1, -^| t=0Q = 0, 



du 

u(t = 00) = 1, -j7|t=oo = 



o"y(t = 00) = 0, — — | t=00 = 
dt 



(32) 



(33) 



With these boundary conditions, integrating equation 
([3"Tjl once over t, we obtain 



d 
di 



5y 



u c dt 2 



= 2 j dq q- 4 Su(q)u c (q) (34) 



B. First-order quantum corrections to the 
Chandrasekhar mass and radius 

In terms of the variable t = 1/x the Chandrasekhar 
mass (fTTf reads 



h = t c ~ tit, T 



Su(t c ) 

I du c \ 
V dt He 



(38) 



Consequently, the dimensionless radius of the gas sphere 



X\ = 



1 1 r\T 



(39) 



is enhanced (diminished) in comparison to the Chan- 
drasekhar radius t~ x when r > (r < 0). 
With the Taylor expansion 



y c (t) ~ y c (t c ) + (^j-) tc (t t c ) + ~(^k(* ~ tcf 
we obtain 



dy c _ ,dy c 



f d 2 y C s 
Vr(^)t c , 



so that the Chandrasekhar mass (1331) becomes 



M(i) = 47Tp c ^ ( (-^) tc +r,B, 



(40) 



(41) 



where 



Af(a?i) =4 7 r / 9 c A 3 — 
at 



?7t 4 d 2 u 
~'dt 2 



(35) 



where t\ = \jx\ denotes the largest zero of the function 
u{t x ) = 0. With the ansatze j26j and ((28j we find to first 
order in 77 



[°° u c {q)5u{q) 6u(t^d 2 y c 

Bc = 2 L dq — (¥k ( ^ k (42) 

The first term in Eq. (|4"Tj) represents the standard Chan- 
drasekhar mass. The quantum correction scales linearly 
proportional to 77. 



5 



C. Approximate solution of the differential 
equation <|3 1 1> 



implying from equation (|48[) 



It is well-known^ that the standard Lane-Emden equa- 
tion, i.e. Eq. with r\ a = 0, is analytically solvable 
only for the values a = 0,1,5 but not for the physically 
interesting case a = 3. However, for general values of a 
the standard Lane-Emden equation has the asymptotic 
solution^ 



Sy(t) 



24i 2 



a 



20t 2 



(52) 



From equation (j2"9")l we obtain accordingly 



4a - 5 



30t 2 



(53) 



y c (t) * i - ^ 



120i 4 



implying 



. , „ in , . a a 



(43) 



(44) 



D. Chandrasekhar mass and radius reductions 

Collecting terms we obtain for the quantum plasma 
correction (j42]) to the standard Chandrasekhar mass the 
quantity 



L 2-(2/ a) „ 1 



a — 1 a(a — 1) 



Gf 2 



l-(2/a) „ j _ 



12i 2 



120i 4 ' 



t(a - 2) 
240i 4 ' 



and 



A 2 



a 
'2*4 



a 

m 2 



respectively. Substituting 



<%(*) = g(t) 



t 4 d 2 U r 



u c dt 2 

yields for the differential equation ((3T 



12i 2 



dt 2+ t* Uc 9 ~ 



,l-(2/q) ^ 2 "c 



rfi 2 



(45) 



(46) 



(47) 



(48) 



(49) 



The boundary conditions f|33|) imply g(t = oo) = a/2. 
Inserting the asymptotics (|4~4l - (l4"Tl) yield for the differ- 
ential equation (j4"9"|) 



dt 2 



f 4 



6i 2 



2 ^ q(a - 1) 
_ ^ 120i 4 



R - a 



5 - 4a 
30t? 



(54) 



Likewise, the corrected normalised radius (13T)1) is 



a 2 n / 5 — a 

1+ — ' 1+ w 



(55) 



For a — 3 the corrected Chandrasekhar mass (|4ip thus 
becomes 



977 
8tf 



1 - 



30< 2 



(56) 



whereas the corrected normalised radius (1551) is 



xi 



tr 



1 



9// 



1 



15£? 



(57) 



Strictly taken, expressions (15^1) - ([57)1 are only valid for 
large values of t c <C 1 because the corrections have been 
derived from the corresponding asymptotic solutions. 
Nevertheless, we use them here fori t c = 1/6.89685 = 
0.145, which is not large compared to unity, in order to 
obtain a first crude estimate of the strength of the quan- 
tum corrections. We find that the quantum corrections 
reduce the standard Chandrasekhar mass by the negligi- 
bly small factor 



a 
2tA 



3a- 2 



I2t 2 



which, at large t, has the asymptotic solution 



(50) 



B c r] = -3.73 • 10 3 ?? = -3.4 • 10" 35 (a 9 x) _1 (58) 

and enhance the normalised stellar radius by the negligi- 
bly small factor 



a a 
2 ~ 480i 4 



(51) 



9rj 



15t? 



= 4.69^ = 4.4- 10- 38 (a 9 x) _1 (59) 
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IV. SUMMARY AND CONCLUSIONS 

The proper quantum plasma treatment of the elec- 
tron gas in degenerate stars such as white dwarfs pro- 
vides an additional quantum contribution to the elec- 
tron pressure. We have calculated how this additional 
pressure term modifies the equation for hydrostatic equi- 
librium, resulting in the quantum modified Lane-Emden 
equation for polytropic equation of states. The addi- 
tional pressure term also modifies the expression for the 
limiting Chandrasekhar mass of white dwarfs. We de- 
velop an approximate solution of the quantum modified 
Lane-Emden equation for general polytropic indices. We 
demonstrate that the quantum corrections reduce the 
standard Chandrasekhar mass by a negligibly small value 
of order 0(1O -35 ) and enhance the stellar white dwarf ra- 
dius by a negligibly small value of order 0(1O -38 ). Devia- 
tions from charge neutrality in the white dwarf electrron- 
ion plasma are of order O(10 -37 ). 

Our study is preliminary in two important aspects. 



First we have used asymptotic expansions of our approx- 
imate solution for the QMLE equation for large values 
of the dimensionless variable t to estimate the quantum 
corrections to the white dwarf mass and radius. Such 
can be improved in future work by a numerical solution 
of the QMLE equation for all values of t. Secondly, we 
have used a nonrelativistic calculation for the additional 
quantum pressure term to combine it with the relativistic 
degeneracy electron pressure term when estimating the 
corrections to the standard Chandrasekhar mass. This 
shortcoming is currently difficult to fix, as a fully rela- 
tivistic kinetic theory of quantum plasmas is not com- 
pletely fleshed out as of yet. 
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dn 
If 



d(nu) 
dz 



= 0, 



(A2) 



and the dynamical electron equation 



du du 
= TT7 + u— 
at dz 



1 d{P c + P Q ) e d$ 



mn dz m dz 

with the Fermi pressure P c and the Bohm pressure 



(A3) 



P 



Q 



2m 



dn 1 ' 2 
dz 



T/2 



d 2 ^' 2 

dz 2 



(A4) 



$ denotes the electrostatic potential and n and rn refer 
to the electron and ion number densities, respectively. 
With the identity 



1 dP Q 

im dz 



h 2 l d 

2m 2 n dz 



dn 1 ' 2 



,1/2 



d 2 n^ 2 
dz 2 



Appendix A: Quantum hydrodynamical equations 

According to Manfredi and Haas& the quantum hy- 
drodynamic equations of an electron-ion plasma in a 1- 
dimensional cartesian geometry are given by the Poisson 
equation 



2m 2 dz 



,1/2 



(A5) 



the dynamical electron equation in full 3-dimensional 
form reads 



dz 



— = Aire [n - m 



(Al) 



the equation of continuity 



du , _ d . _ e ^ ^ 1 _ _ 

-57 + (it - 7^)u+ — V$ + VP C 

dt dx m mn 



h 2 /V^i/2 
2m 2 { nV2 



(A6) 



7 



1. Quantum modified hydrostatic equilibria 

For static (u — 0) systems, after adding gravity, the 
electron dynamical equation reads 



= mng(x) + enV$ + VP C - 




Because the ion contribution to the pressure is negligibly 
small, the corresponding ion equation reads 

~ mi-nigix) - Zen t V<b (A8) 

Adding equations (|A7|) and (| A8|) and using the charge 
neutrality condition Zrii — n eliminates the electrostatic 
potential leaving 



where we have inserted Eq. (|A9|) . This electric field pro- 
vides a net space charge density a that can be estimated 
by calculating the divergence of Eq. (|A12|) as 

-* TYl ' 

V • E = 4na = Aite(Zn t - n) = -V • g (A13) 

eZ 

With V -g = —AttGp and p = we obtain the positive 
space charge density 

a=e(Zm-n) = ±-i = AH 

eZ eZ z 

that has to be screened by an appropriately negative 
space charge density on the surface of the star—. As rel- 
ative deviation from the overall charge neutrality state 
we obtain 



. (run + + Vfi - ^ j (A9) |^ _ ^ _ £ft _ 1.01 . 10"" (A15) 

With the mass density p — m^n, + mn ~ rriin/Z we if we use Z = 2 and m; = m p . The deviation from the 
derive overall charge neutrality state is negligibly small. 



V 2 n x / 2 \ Zm-n Zm-n Gm 2 



^ i h 2 z /v 2 p 1/2 

= g(x) + -VP C - V 

p Zmrrii \ p L / z 



(A10) 



which is identical to Eq. (|4|). Under charge neutrality 
the electric field E = — V<i> does not enter the equation 
for stellar hydrostatic equilibrium. 



2. Electric field equation and charge neutrality 
state 



Multiplying equation (|A8[) with Z and subtracting 
equation (|A7[) yields for the electrostatic potential^ 

(Z^erii + en) V<£> = (Zm^i — mn) g(x) 



h 2 n_ /V 2 nV2 



2m 



,1/2 



(All) 



Because of the stronger action of the gravitational field 
on the plasma ions as compared to the plasma electrons, 
a charge separtion results that induces a nonzero electro- 
static potential and electric field. Assuming only small 
perturbations of the charge neutrality condition Zni ~ n 
then yields 



E 



n 2 



e(Z + 1) 

W/2 



e(Z + l)n 



2me(Z + 1) 



,1/2 



= -^(f) (A12) 
eZ 



Appendix B: Relativistic generalisation of the 
quantum modified Lane-Emden equation 

Assuming only isotropic stresses, the mixed compo- 
nents of the energy-momentum tensor of a spherically 
symmetric distribution of fluid obeying a polytropic 
equation of state in a co-moving system read 



T, 1 = T 2 = T 3 3 



-P 



n 



pc 



(Bl) 



where P is the pressure and e = pc 2 is the energy density. 
Using a spherical form of the metric at rest with respect 
to the fluid distribution 



ds 2 = e^c 2 dt 2 - r 2 (d6 2 + sin 2 6d<t> 2 ) - e x ^dr 2 (B2) 
we obtain the time-independent gravitational equations 7 



\dv 1 \ 1 8ttG 
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8ttG o 




e 
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dr 


r 2 ) 


+ r 2 = C 4 PC 


(d 2 v 
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dv 


2 1 


dv 


d\~ 


1 dv d\\ 


\dr 2 ' 


f 2 


dr 


+ - 
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dr 


dr 


~2~dr~dr) 



(B3) 



(B4) 



16ttG 
(B5) 
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Combining these equations we derive the expression 



order C(ft 4 /,4 4 ) 



dP 
dr 



P- 



pc 2 ) 



dv 
dr 







(B6) 



which is the general relativistic analogue of the classical 
form of equation ((5]). Introducing the quantum correc- 
tion Pq in the expression for the pressure P = P c + Pq 
(see equation (A18) and (A20)) equation (B6) yields 



dw l-2tt+lw/( 2 4 
w + vy a (-7T + — C -7T 

«C 1 + vy a rfC 



v a J 2 ^\v a J 2 



dya 



I 2 \ 2 



d( 



1 dP r 



d /V2pV2( r ) 



p(r) dr Imrrii dr \ p x l 2 (r) 
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(B8) 
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Equation (|B4|) yields with the definitions and substitu- 
tions given previously 



Using the polytropic equation of state (|12p and the 



• C, where A is de- 



substitutions w(Q = —, — 

2(q + 1)<7 

fined by expression (fT5|) . C = ^/-A and the constants 

h 2 

Xa 



2n»MV(l + ' " = ^ /7c2 ^ inSerting 
equation (|B7|) into (|B3|) we derive neglecting terms of 



— = CV 



(B9) 



In the non-relativistic limit the general relativistic quan- 
tum system (|B8|) and (|B9|) reduces to equation ([T4|) , while 
for /i — > it reduces to the general relativistic Lane- 
Emden equation of stellar structure^. 



